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We investigate the ground state of interacting spin-| fermions (3D) at a finite density (p ~ k%) in 
the presence of a uniform non-AbeUan gauge field. The gauge field configuration (GFC) described 
by a vector A = (A^;, Ay, Az), whose magnitude A determines the gauge coupling strength, generates 
a generalized Rashba spin-orbit interaction. For a weak attractive interaction in the singlet channel 
described by a small negative scattering length {kF\as\ < 1), the ground state in the absence of the 
gauge field (A = 0) is a BCS (Bardeen-Cooper-Schriefi^er) superfluid with large overlapping pairs. 
With increasing gauge coupling strength, a non-Abelian gauge field engenders a crossover of this 
BCS ground state to a BEC (Bose-Einstein condensate) ground state of bosons even with a weak 
attractive interaction that fails to produce a two-body bound state in free vacuum. For large gauge 
couplings (X/kp S> 1), the BEC attained is a condensate of bosons whose properties are solely 
determined by the gauge field (and not by the scattering length so long as it is non-zero) - we 
call these bosons "rashbons". In the absence of interactions (as = 0~), the shape of the Fermi 
surface of the system undergoes a topological transition at a critical gauge coupling At. For high 
symmetry gauge field configurations we show that the crossover from the BCS superfluid to the 
rashbon BEC occurs in the regime of A near At- In the context of cold atomic systems, this work 
makes an interesting suggestion of obtaining BCS-BEC crossover through a route other than tuning 
the interaction between the fermions. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.-d, 67.85. Lm, 71.70.Ej 



I. INTRODUCTION 

Recent experimental progress in the generation of syn- 
thetic gauge fieldJl^ has enhanced the possibilities of 
controlled experimental studies of outstanding problems 
of quantum condensed matter and even of high energy 
physic^ using cold atomic systems. Many theoretical 
works^ * have explored the possibilities of generating 
both Abclian and non-Abeilan gauge fields. The experi- 
mental work with synthetic gauge fields have been with 
bosonic ^''Rb atoms. In a recent commentary -' , the in- 
vestigation of fermions in synthetic non-Abelian gauge 
fields has been identified as a key research direction. 

The study of interacting fermions in 3D space with 
controlled interactions has been one of the key suc- 
cesses of cold atoms research.l^^HUl A particular example 
is the problem of the crossover from a BCS ground state 
with large overlapping pairs to a BEC of tightly bound 
bosonic pairs with increasing strength of attractive inter- 
actions - a p henomenon that was suggested many years 
earlier .n^HlS] The superfiuid transition temperature on the 
BCS side is determined by the superfluid energy gap, 
while on the BEC side the transition temperature is de- 
termined by the condensation temperature of the tightly 
bound bosonic pairs of fermions.'^ A review of BCS-BEC 
crossover that is particularly useful in the context of this 
paper may be found in reference |17| . 

An interesting question is regarding the fate of inter- 
acting fermions in the presence of a non-Abelian gauge 
field. Motivated by the fact that even a spatially uniform 
non-Abelian ga uge field produces interesting physical ef- 
fects for bosons ,121121121 -^^e focus on interacting fermions 



in uniform non-Abelian gauge fields. 

In a recent paper,'2Sl two of us investigated how a uni- 
form non-Abelian gauge field influences the bound state 
of two spin-^ fermions interacting via a contact attrac- 
tion in the singlet channel characterized by a s-wave 
scattering length a^. The type of uniform non-Abelian 
gauge field considered in that work leads to a general- 
ized Rashba spin orbit interaction. A key finding of that 
work is that for high symmetry gauge field configura- 
tions (more precisely defined in the next section), a two- 
body bound state exists for any scattering length however 
small and negative. The study suggested that the BCS- 
BEC crossover is drastically affected by the presence of 
a non-Abelian gauge field. 

Here we pursue the ideas of reference [2^ through a 
study of the ground state of a finite density of interact- 
ing fermions in a non-Abelian gauge field by means of 
mean field theory. We show that increasing the strength 
of a non-Abelian gauge field produces a crossover from 
a BCS superfiuid (which is the ground state in the ab- 
sence of the gauge field) to a BEC of bosons at a fixed 
interaction (fixed scattering length a^) however small and 
negative. Further, the bosons that condense to form the 
BEC at large gauge couplings are tightly bound pairs of 
fermions whose properties are determined solely by the 
non-Abelian gauge field ~ we have ventured to call these 
bosons "rashbons" . For a given attractive interaction 
(fixed Og), therefore, the crossover takes place the stan- 
dard BCS superfluid state to a rashbon BEC. There is 
an additional feature of the crossover that is particularly 
noteworthy. The Fermi surface of the non-interacting 
system (os = 0~) undergoes a transition in its topol- 
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ogy with increasing gauge coupling strength. We show 
that for high symmetry gauge field configurations, the 
crossover regime of gauge couplings in the presence of 
interactions overlaps with the regime of topological tran- 
sition of the non-interacting Fermi surface. In a sense 
this provides a "geometrical" view of the crossover. 

The statement of the problem we address and a de- 
tailed summary of our results are given in Section [iT] 
The mean field formulation is detailed in Section HTll Sec- 
tion |IV| describes the results. The paper is concluded 
with a discussion in Section |Vl We recommend the read- 
ing of Section [ll] and Section |V] to obtain a physical pic- 
ture of our results. 



II. QUESTION ADDRESSED AND SUMMARY 
OF RESULTS 

In units where the mass of the fermions and Planck's 
constant are set to unity, the hamiltonian of the fermions 
moving in a uniform non-Abelian gauge field is 



-Rgf = j dV«'^(r) 



(pi - A'^T^) • (pi - A^T^ 



vl'(r), 
(1) 

where ^'(r) = {tpfj{r)},(j =t,i are fermion operators, p 
is the momentum, = A^Bi, arc uniform gauge fields, 
(/i = x,y,z) are Pauli matrices and e^'s are the unit 
vectors in the i-th direction, i = x,y,z. As in |20j . we 
specialize to ~ \i6^ leading to a hamiltonian with a 
generalized Rashba spin-orbit interaction 



(2) 



where px = J^iPi^i^i- The vector A = AA = J2i^i^i 
describes a gauge field configuration (GFC) space as de- 
picted in fig. [T] we call A = |A| as the gauge coupling 
strength. High symmetry GFCs are important and have 
been classified in [5D] as prolate, spherical and oblate. Of 
particular interest are the configurations shown in fig. [T] 
called extreme prolate (EP), spherical(S) and extreme 
oblate(EO). 

The one-particle states of Hr are 



that disperse as 



\ka) = |fe) <g) \akx) 



£fea = ^ - a\kx 



(3) 



(4) 



where fc-the momentum, and a = ±l-the eigenvalues 
of the helicity operator px ■ t, are the good quantum 
numbers. The quantity kx is defined analogously with 
Px m cqn. (|2|. For any A, the two helicity states for a 
given k are degenerate only at fc = 0. 

The interaction between the fermions is described by 
contact attraction in the singlet channel 



n^^v d^r^\{r)i;lir)i;^{r)ij^{r). 



(5) 




FIG. 1. (Color online) Gauge field configuration (GFC) space. 
The non-Abelian gauge field of eqn. ([2| is described by a 
vector A = {X^, Xy, X^) = AA, where A — |A| is the gauge 
coupling strength and A is a unit vector. High symmetry 
GFCs such as extreme prolate (EP, A — (0,0,1)), spherical 
(S, A = ^(1, 1, 1)) and extreme oblate (EO, A = ^(1, 1,0)) 
are as shown. 



The endemic ultraviolet divergence of the theory de- 
scribed by the hamiltonian 



(6) 



is handlecP^ by introducing an ultraviolet momentum 
cutoff A. This entails characterization of the attraction 
by a physical parameter that describes the low energy 
scattering properties while making the parameter v de- 
pend on A. More precisely. 



1 



+ A 



1 

47ra, 



(7) 



where is the s-wave scattering length in free vacuum, 
i. e., when the gauge field is absent (A = 0). In free 
vacuum (3D) only an attraction larger than a critical 
strength can produce a two-particle bound state. This is 
embodied in the fact that for Os < (BCS side) there 
is no two-body bound state; a bound state develops only 

as Qs —i' — oo, or > 0~ (resonance). For > (BEG 

side) a bound state is obtained with a binding energy 
Eh — ^ • 

A uniform non-Abelian gauge field brings about re- 
markable changes in the two body problem as shown in 
reference [20]. Most vividly, for high symmetry GFGs 
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TABLE I. Summary of the two body proble nPin high symmetry GFCs (A > 0). Et is the binding energy and r]t is the 
triplet fraction (not reported in reference [2^). The bi-axial spin nematic structure is similar to the BW (Balian-Werthamer) 
or B-phase of "^He, and the uni-axial nematic structure to that of the ABM (Anderson-Brinkman-Morel) or A-phase of '^He. 
The values of Et and rjt at resonance, which correspond to the properties of the rashbon (for S and EO cases), are exact results, 
while others are asymptotic values. 



such as S and EO, the critical scattering length a^j, re- 
quired for the formation of a bound state vanishes, i. e., 
there is a two body bound state for any scattering length 
however small and negative (deep BCS side). The size 
of the binding energy of the bound state depends on the 
GFC. On the BCS side, the binding energy has an ex- 
ponential dependence on the scattering length and A for 
the EO GFC while for the S GFC this dependence is 
algebraic. Another interesting aspect that emerges is 
the symmetry of the bound-state wave function. In a 
non-Abelian gauge field the normalized bound-state wave 
function is made up of spatially symmetric singlet and 
spatially antisymmetric triplet pieces 

\A) = m + \^t). (8) 

Time reversal symmetry of the hamiltonian is preserved 
and this two-body-bound-state wave function picks up 
a nematic spin structure consistent with the symmetry 
of the GFC .122 The triplet content of the wave function 
which also measures the "amount of nematicity" is char- 
acterized by a parameter which we call the triplet fraction 

Vt = i^tlA)- (9) 

The binding energy £';,, the triplet fraction ijt and the 
spin symmetry of the two-body wave function for differ- 
ent GFCs are summarizecP^ in Table |lj The gist of refer- 
ence is that high symmetry GFCs induce high degen- 
eracy in the low energy (infrared) one particle density of 



states and this promotes bound state formation. Collo- 
quially, high symmetry GFCs are "attractive interaction 
amplifiers" . An aspect of the two body problem that is 
important in the discussion below is that the physics of 
the two body problem in the presence of the gauge field 
(A > 0) is completely determined by the dimensionless 
parameter Acs . All aspects of the solution depends only 
on Xag when length and energy are respectively measured 
in units of A""'^ and A^. This is true for any GFC except 
EP GFC as is evident from Table. U 

In this paper we investigate the system described by 
the hamiltonian of eqn. ([6| at finite density p of the 
fermions. The finite density of particles introduces an 
additional energy scale which can be conveniently taken 
to be the Fermi energy Ep (and an associated Fermi wave 
vector kp) in the absence of the gauge field (A = 0) 

Ep='f = l{3.'pf\ (10) 

At finite densities of fermions, therefore, the ground state 
of the system (zero temperature) of eqn. ([6| is determined 
by the dimensionless parameters kpag, X/kp, and the 
direction A in GFC space. In this paper we work at fixed 
density and study the evolution of the ground state of the 
system with X/kp and kpag for various As corresponding 
to high symmetry GFCs. 

The possibility of interesting physics in this system is 
suggested by the following observation which provides 
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FIG. 2. (Color online) Fermi surface topology transition (FSTT) with increasing gauge coupling strength for EO GFC. (a) 
Two overlapping spherical Fermi surfaces in the absence of a gauge field (b) A gauge coupling smaller than At: The union of 
the + and — Fermi surfaces forms a spindle torus. The apple of the spindle torus is the + helicity Fermi surface which is shown 
with a blue border, and the lemon of the spindle torus is the — helicity Fermi surface which is shown with a red border, (c) 
The gauge coupling that obtains the FSTT. The + helicity Fermi surface is a horn torus, while the — helicity Fermi surface 
vanishes, (d) For a gauge coupling larger than At, there is only the + helicity Fermi surface which is a ring torus. Note that 
all the figures show only a sectioned half of the Fermi surfaces. 



the motivation for this work. Consider a system of non- 
interacting (NI) fermions (T-L^ — in eqn. Q). In the 
absence of a gauge field (A = 0), the ground state has a 
chemical potential Ep and is described by two identical 
filled Fermi seas bounded by spherical Fermi surfaces of 
radius kp - one each for f and J, spins. In the presence 
of the gauge field (A ^ 0), the helicity a is the good 
quantum number along with momentum (see eqn. (|3|), 
and hence the ground state will be two Fermi seas, one 
for each helicity. The chemical potential-- now depends 
on A through a function fi^iW that is determined by 
A. Both of the Fermi seas are generically non-spherical 
with a shape determined by A. Since the one particle 
states with opposite helicities but with same momentum 
are non-degenerate (see eqn. ([4])), the Fermi surfaces of 
different helicities are not identical and evolve differently 
with increasing A (at a fixed density p). The most inter- 
esting aspect is that since the -I- helicity state is lower 
in energy than the — helicity one for all momenta, upon 
increasing A, the volume enclosed by the -I- helicity Fermi 
surface increases at the expense of that of the — helic- 
ity Fermi surface. Matters come to a head at a critical 
gauge coupling A^ (which depends on A) where the — 
helicity Fermi sea ceases to exist since the chemical po- 
tential fJ.NiW falls below the bottom of the — helicity 
band. Thus, for A > At the ground state is a Fermi sea 
of only -I- helicity. This is illustrated for the EO GFC 
in fig. [2] The values of At determined by the density 
of particles for various high symmetry GFCs are given 
in Table |llj it is to be noted that in all cases At is of 
order kp. Another aspect to be noted is that there is a 
change in the topology of the -|- helicity Fermi surface at 
At. We call this the Fermi surface topology transition 
(FSTT) - hence the subscript T in At- For example, in 
the EO case, the genus of the -I- helicity Fermi surface 
changes from zero (homeomorphic to a sphere) to unity 



(homeomorphic to a torus) at At as illustrated in fig. [2j 

What happens when the interaction (eqn. ([5|) is 
turned on? Consider an interaction with a small nega- 
tive scattering length {kp\as\ <C 0, deep BCS side). For 
A <C At, the ground state is a superfluid with overlapping 
pairs with an exponentially small excitation gap. The 
chemical potential of this state is nearly unaffected and 
is /ijv/(A), that of the non-interacting system in a gauge 
field. The only qualitative difference from the usual s- 
wave BCS state is that the pair wave function now has 
a small triplet content and associated spin nematicity 
induced by the gauge field. This picture changes dras- 
tically in the case of high symmetry GFCs (such as S 
and EO) when the gauge coupling strength A is tuned 
past At- The key finding of this paper is that for high- 
symmetry GFCs, a BEG oj tightly hound pairs is obtained 
for A 3> At even with a small negative scattering length. 
In other words, one can engineer a BCS-BEC crossover 
with a high-symmetry GFC by increasing the gauge cou- 
pling strength even with a very weak attractive interac- 
tion that is unable to produce a two-body bound-state 
in free vacuum. This result arises from the fact that 
for A 3> At, the size of the two-body bound-state wave 
function (see Table |l| becomes smaller than the inter- 
particle spacing. The fermions therefore form tightly 
bound pairs which then Bosc condense in the zero center 
of mass momentum state. As is evident from the discus- 
sion, the physics of these results owes to the character 
of high-symmetry GFCs to act as attractive-interaction- 
amplifiers. Indeed, as A/At — > oo the chemical potential 
tends to that determined by the two-particle bound-state 
energy. Since A|as| — > go (fixed as), the nature of the two 
body hound state obtained is identical to that obtained 
with a resonant scattering length in the presence of the 
gauge field as tabulated in Table [Tj The properties of this 
bosonic bound state of two fermions is determined solely 
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TABLE II. Summary of the properties of the superfluid ground state for two regimes of the gauge coupling strength A. The 
results shown here are for a weak attractive interaction (as < 0, fcF|as| ^ 1). is the ground state chemical potential and rjt 
is the triplet content of the pair wave function. Well before the FSTT (A <^ At), the ground state is a BCS superfluid. In the 
regime (A S> At ), for the S and EO GFCs, the chemical potential is close to half of the rashbon energy (see Table. |I| and the 
pair wave function attains the same triplet content as that of the rashbon, clearly indicating a crossover from a BCS superfluid 
to a rashbon condensate. Any CFC except EP produces a BCS - rashbon BEC crossover. 



hy the Rashha gauge field] we call these emergent bosons 
as "rashbons" (see Sec. IV B and second para of Sec. |V] 
for details). The BEC that is obtained for A ^ At is 
a rashbon condensate. The results for various GFCs are 
tabulated in Table [n] which is a summary of this paper. 

In the remaining sections we illustrate these conclu- 
sions by a mean field theory of the superfluid ground 
state of this interacting fermion system. Mean field the- 
ory is known to give a qualitatively correct description 
for the superfluid ground state.l^ 



III. MEAN FIELD THEORY 

We now describe the details of the mean field anal- 
ysis of the superfluid ground state of fermions in a 
non-Abelian gauge field. This analysis involves cer- 
tain straightforward manipulations beyond the standard 
formulation,^^ and hence presented in detail. We note 
that the present analysis can treat any GFC. Results for 
specific GFCs of interest will be presented in the next 
section. 

To perform a mean-field analysis of the superfluid 
ground state we recast the interaction term 'H^J in a con- 
venient form 



= - / d?rS\r)S{r), 



(11) 



where {r) is the singlet creation operator 

SHr) = [^{r)i,l{r) - V4(r)^|(r)) . (12) 
The interaction in terms of the singlet operators can be 



cast in momentum space 



(13) 



where V is the volume of the system and the mean-field 
ansatz corresponds to taking 



{S{q)) = {S)S,^o, 
where (S) = (S'(O)) with 



5t(0) = ^V 



(14) 



(15) 



The fermion operators c],^ are defined by 



-'fccr 



d^re-'''-''tpl{r). 



(16) 



It is now convenient to sum only over half of the allowed 



k values in eqn. ( 15 ) 



5(0) = V2^ 



(17) 



as is here and henceforth indicated by the prime over 
the summation symbol. The advantage of this exercise 
is that the operator 5(0) can be written in the helicity 
basis as 



(18) 
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Introducing a chemical potential fj,, we obtain the mean- 
field Hamiltonian as 



Hmf = CfeacLcfca + A ^ «4aclfca 

ka koL 



where cj[.^ are electron operators associated with the one- 
particle helicity eigenstate (eqn. (3|), A = ^^jS^ is the 
order parameter (taken to be realj7and £^ka = £ka — A*- 
Here is £ka referred to the bottom of the -I- helicity 



band. Noting inversion symmetry, S,-ka — £,kai eqn. (19) 
can now be recast as 



ka 



^kc 



-koL 



(20) 



ka 



which now has the standard form except for the fact that 
the summation over k is carried out only over half of the 
momentum space and a sum over the two helicities is 
taken. 



The hamiltonian in eqn. ( 20 1 can now be diagonalized 
in terms of the Bogoliubov quasiparticle operators as 

■Hmf = J2 ^'"^ ('^Li7fcai + 7L2 



ka 



(21) 



ka 



where Eka — \/ + A^ (A is also the excitation gap), 
and 



7fcal = UkaCka - OiVkc[_^^ 
7L2 = O^VkaCka + UkC^-ka 



with 



^ka 



; 1 



^ka 
Eka 



^ka^ ^[^~ -n^ 
^ \ rjka 



(22) 



(23) 



A standard analysis now leads to the gap equation 

-- = y— . (24) 



ka 



Noting the inversion symmetry of the problem and using 
the renormalization of the interaction, the gap equation 
becomes 



ka ^ 



(25) 



The number equation is 



-E- 



V 



^ka 
Eha 



(26) 



(19) The solution of eqn. (|25|) along with the number equation 



eqn. (26), determines the chemical potential and the 



gap parameter A in the ground state. 

The ground state |^g) of the system is given by 

I*g) - n^"-^" + J|0) (27) 

ka 

where |0) is the fermion vacuum. This can be (up to a 
normalization) be re-written as 

|*G)=e^>) (28) 
where is the pair creation operator given by 



-ka 



(29) 



where (jyua 



^ka 



The singlet and triplet parts of the 



pair can be extracted by noting that 



ka 



singlet 



(30) 



triplet 



with 



,(fe) = 
tik) = 



(31) 



This analysis sheds light on how an attraction in the 
singlet channel in presence of a non-Abelian gauge field 
can produce a triplet piece in the pair wave function. 
The triplet content rjt is now defined as the weight of the 
triplet piece of the pair creation operator in eqn. (30 1. 



One can also charecterize this by an expectation value 
of the quadrupole operator of reference [22J. However, 
this definition for the triplet content rjt given above pro- 
vides a physically transparent and a simple measure of 
the quantity of interest. 

A remark about the Bogoliubov quasiparticles ob- 
tained in eqn. (22) is in order. It appears that for each 



helicity there are two branches of quasi-particle excita- 
tions labeled 1 and 2. Ostensibly, therefore, there are 
four branches of quasiparticles which at the first sight is 
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FIG. 3. (Color online) Evolution of the triplet content r/t of 
the pair wave function as a function of the gauge coupling 
strength A for an EP GFC with kpUs = —1. The evolution of 
the same quantity of the non interacting system {as = 0~) is 
also shown for comparison. 



surprising. Note, however, that these four branches are 
defined only in half of the momentum space. If the Bo- 
goliubov excitation were defined for all fc, they will not 



t -^t 
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be independent, for example, 
motivation behind introduction of the sum over one half 
of the momentum space in eqn. (17). It is now clear that 



the formulation recovers the correct counting of excita- 
tion, i.e., within the present formulation, two excitations 
for each k in momentum space is recovered as four exci- 
tations for each k in half the momentum space. 
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IV. RESULTS FOR SPECIFIC GAUGE FIELD 
CONFIGURATIONS 



In this section, we shall present results of how the 
ground state of the system evolves with A for various high 
symmetry GFCs. We shall be concerned only with neg- 
ative scattering lengths (a^ < 0) since this is the regime 
which has the most interesting physics. In the absence of 
the gauge field (A = 0) there is no two body bound state, 
and for as < the usual BCS superfluid ground state^'' 
\BCSq) is obtained. For small A, i. e., A <C At, we expect 
and find the ground state I^'q) to be qualitatively close 
to \BCSo) state with an exponentially small excitation 
gap and a chemical potential essentially unaltered from 
that of the non-interacting problem iini{X). When A is 
increased beyond Ay, we find, in some cases (S and EO 
GFCs), that the chemical potential begins to fall and 
approaches —Eb/2, the value set by energy of the two 
body bound state. This signals the crossover to the EEC 
state. Additionally, the pair wave function defined by 



eqn. ( 29 ) approaches the wave function of the two body 



bound state. 

A summary of the results for various GFCs discussed 
below is given in Table. [ll| 



FIG. 4. (Color online) Evolution of the ground state of a 
collection of interacting fermions {kpas = — j) with gauge 
coupling strength A for the S GFC. (a) Chemical potential ob- 
tained from a numerical solution of mean field theory (MFT) 
is compared with the chemical potential of the non interact- 
ing system (NI) and that set by the binding energy of the two 
body problem (— i?i,/2). For A < At the chemical potential 
is indistinguishable from the that of the non-interacting sys- 
tem. For A > At the chemical potential approaches the two 
body value indicating a crossover to a BEC. (b) Evolution of 
the numerically obtained mean field energy gap A with the 
gauge coupling strength A. The analytical result eqn. (361 



is also shown and is indistinguishable from the numerical re- 
sult, (c) The dependence of the triplet content (774) of the 
pair wave function defined in eqn. ( 29 1 on the gauge coupling 
strength. This is compared with the same quantity of the 
non-interacting system (NI) and with that of the wave func- 
tion of the two-body bound state. It is seen that the pair wave 
function evolves to two-body bound-state wave function. 



A. Extreme prolate (EP) GFC 

This GFC with A = (0, 0, A) has an FSTT at At = kp- 
Before FSTT (A < At), the -I- helicity Fermi sea consists 
of the volume enclosed by two intersecting spheres of ra- 
dius kp centered around (0,0, ±A), while the — helicity 
Fermi sea is the lens shaped region formed by the volume 
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common to both spheres. When A exceeds At the — he- 
Ucity Fermi surface vanishes, and the + hehcity Fermi 
sea is made of two disjoint spheres centered at (0, 0, ±A). 
The chemical potential //at/ (A) = Ep, i. e., is unaffected 
by the EP gauge field. 



EO, kpa 




FIG. 5. (Color online) Evolution of the ground state of a 
collection of interacting fermions {kpCLa = —1) with gauge 
coupling strength A for the EO GFC. (a) Chemical poten- 
tial obtained from a numerical solution of mean field theory 
(MFT) is compared with the chemical potential of the non 
interacting system (NT) and that set by the binding energy 
of the two body problem {—Eh/2). For A < At the chemical 
potential has qualitative features similar to those the non- 
interacting system - the numerical differences are due to the 
larger value of the magnitude of the scattering length. For 
A > At the chemical potential approaches the two body value 
indicating a crossover to a BEC. (b) Evolution of the numer- 
ically obtained mean field energy gap A with the gauge cou- 
pling strength A. The analytical result closely approximates 
the numerical result, (c) The dependence of the triplet con- 
tent {rit) of the pair wave function defined in eqn. (291 on 
the gauge coupling strength. This is compared with the same 
quantity of the non-interacting system (NT) and with that of 
the wave function of the two-body bound state. It is seen that 
the pair wave function evolves to two-body bound-state wave 
function. 



For feplflsl <^ 1, the standard resullP^for the excitation 
gap is 

A 8 ^ 



and the chemical potential is 

fifv Ef 



(32) 



(33) 



Not unexpectedly, the excitation gap A and the chemical 
potential are unaltered with increasing A. The ground 
state for any A is a superfluid state with large overlapping 
pairs, and there is no BCS-BEC crossover for the EP 
GFC. There is, however, a qualitative change in the spin 
structure of the pair wave function. With increasing A, 
the pair wave function develops a triplet content r]t (see 
hg. 3l which attains a value close to | at A = At and 
stays so with further increase of A. 

The physics behind this result can be traced to the fact 
that for the EP GFC the kinetic energy content inside the 
non-interacting Fermi sea is unaltered by the increase of 
A. Therefore, the gauge coupling A stays neutral in the 
competition between kinetic energy and the attractive 
interaction. This, again, is the reason why the energetics 
of the two-body problem is unaffected by the presence of 
an EP gauge field (see Table. |l]). This is a feature specific 
only to EP GFCs. 

It must be noted that the non-interacting ground state 
also has a triplet content (see Fig. [s]) . As is evident (see 
eqn. (30)), this arises from the fact that the -I- helic- 
ity Fermi sea is different (and larger) than the — helicity 
Fermi sea. The triplet content of the non-interacting sys- 
tem increases monotonically with A and attains a value 
of i at A = At and remains at this value for any larger A. 
As expected, in the presence of an attractive interaction 
in the singlet channel (os < 0), the pairs have a triplet 
content less than that of the non interacting system. 

We note that the qualitative nature of the results for 
negative scattering lengths (a^ < 0) of larger magnitude 
are similar to those for feFl^sl ^ 1- 



B. Spherical (S) GFC 

When A = ^(1, 1, 1) a spherical (S) GFC is obtained. 
Starting from two identical overlapping spheres at A = 0, 
the non-interacting Fermi surfaces of the two helicities 
continue to be spheres with their centers at the origin 
of the momentum space for < A < At. Here At = 
\/3 

-^^kp- When A ^ At, the chemical potential of the 
non-interacting system depends on A as 



A<jv/(A) 
Ef 



1 

2l 



A 

At 



(A < At) (34) 



In this regime, the radius of the -f helicity Fermi surface 
is larger than that of the — helicity Fermi surface. At the 
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FSTT, the — helicity Fermi surface vanishes and ceases to 
exist for all A > At- After the FSTT, the + hehcity Fermi 
sea is "a sphere with a hole" , i. e., the region bounded by 
two concentric spherical Fermi surfaces. For A ^ At the 
chemical potential of the non-interacting system goes as 



Mjvj(A) 



2S 



T 



(A > At). 



(35) 



Consider now the situation when fci?|as| <^ 1. When 
A = 0, the usual BCS state with properties given by 
eqn. (32 1 and eqn. (33) is the ground state. For A ^ At, 
fj, is very nearly equal to that given by eqn. (34); the 
gap equation can be solved analytically in this regime to 
obtain 



A = 



8^Ar/(A) 



exp 



(A)+A2 



exp 



y2|a,|(6/ijv7(A) + A2) 



(36) 

Fig. |4][a) and (b) show, respectively, the numerical so- 
lutions of the chemical potential and gap as a function of 
A. Fig. |4|^a) also shows the non-interacting chemical po- 
tential, and the two-body energy —Ef)/2 (which depends 
on A and as only). As is evident the chemical poten- 
tial fi is identical to the non-interacting value uniW for 
A ^ At. There is also excellent agreement for the gaps 
obtained from the numerical solution with the analytical 
result given in eqn. ( 36 ) . When A reaches At the chemical 
potential begins fall below fiNi, and on further increase 
of X (A > At ), the chemical potential tends to that set 
by the two body problem. This clearly signals a crossover 
from the BCS like state for A ^ At to a EEC state where 
the fermions from tightly bound bosonic pairs which then 
condense in the zero center of mass momentum state. 

Further corroboration of the crossover to the BEC like 
state with increasing A can be obtained by a study of the 
triplet fraction rjt which is shown in Fig. |4] Again, rjt cor- 
responding to the non-interacting system monotonically 
increases and attains a value of ^ at At. The triplet con- 
tent of the superfluid pair, as expected, is less than that 
of the non-interacting system, but has a similar quali- 
tative behavior as the NI case in the regime A ^ At. 
The triplet fraction attains a maximum at a A close to 
At and then begins to fall. On further increase of A, rjt 
approaches that of the two-body bound-state wave func- 
tion, demonstrating again that the pair wave function 
tends to the two-body bound-state wave function. We 
also see that A = At marks the crossover point, i. e., 
the crossover regime is precisely the regime of A where 
change in the topology of the non-interacting Fermi sea 
takes place. 

It is particularly interesting to study the BEC state 
that is attained when A — ?> c». The key point as noted in 
section, [njis that the physics of the two-body bound state 
is determined by the dimensionless parameter Aa^ (see 

Table. U^. Therefore, as A — >■ oo, the parameter i^^^^ >■ 0. 

Thus tne state that is obtained is same as that obtained 
for the two-body bound state with a resonant scattering 



length in the presence of the gauge field (\ > 0) (Ta- 
ble. |l])! Therefore the properties of the BEC for A — > oo 
are completely determined by X, independent of the scat- 
tering length (as long as it is non vanishing), i. e., the 
system is a collection of Bosons whose properties are de- 
termined solely by the Rashba interaction. Hence we 
call this tightly bound bosonic state of two fermions as 
"rashbon" . Rashbon is a bound state of two fermions 
in a Rashba gauge field (A > 0) at resonant scattering 
length (^ = 0). 

Again, for scattering lengths of larger magnitude, the 
qualitative physics remains identical. We shall illustrate 
this point in the next section by considering the EO case 
with a scattering length of larger magnitude. 



C. Extreme oblate (EO) GFC 

The evolution of the non interacting Fermi surfaces for 
this GFC (A = ^(1, 1,0) is shown in fig. [2] The non- 
interacting chemical potential in the regime A ^ At is 



IJ-NiW 
Ef 



= 1 - 



A 



(A < At) (37) 



and that in the regime A ^ At is 
fiNiW M A ^ At 



Ef 



3tt 



A 



(A > At). 



(38) 



In the regime when A ^ At and for fci?|as| ^ 1, 
the chemical potential is well approximated by the non- 
interacting value. Further, we can obtain a rather lengthy 
analytical expression for the gap (not shown). 

To illustrate that the qualitative nature of the transi- 
tion is unaltered by the size of the scattering length, we 
study this GFC with fcTla^l = 1. The results are shown 
in fig. [5] These results clearly illustrate a crossover from 
the BCS like state to a BEC state of rashbons. Note, in 
particular, that rjt of the many body pair wave function 
tends to that of the two body bound state wave func- 
tion with a resonant scattering length (rashbon) given in 
Table. U 



V. DISCUSSION 

We conclude the paper with further discussion of our 
results. On the BCS side fci?|as| <C 1 and A <C At, the 
transition temperature will be determined by the zero 
temperature gap which we have calculated in this paper. 
On the rashbon BEC side, the transition temperature will 
be determined by the mass of these emergent boson^^ 
which will be renormalized from the value of twice the 
fermion mass due to the gauge field. 

As noted earlier, the rashbon is a bound state of two 
fermions in a Rashba gauge field (A > 0) when the s- 
wave scattering length is infinity, i.e., at resonance. This 
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two-fermion bound state exists for all GFCs except the 
EP GFC and has a spin structure determined by A of 
the GFC (see the "resonance" column of Table. |l]). As is 
evident this state is not rotationally symmetric - it is an 
"anisotropic particle" that emerges. It is also interest- 
ing to contrast the rashbon state obtained in a Rashba 
gauge field (A > 0) with the two-body quasi-bound state 
obtained in free vacuum (A = 0) at resonance. In the lat- 
ter case, the binding energy is zero, and the state is scale 
free with a singlet spin structure. This is to be contrasted 
with the rashbon state whose binding energy is A^ times 
a dimensionless number that depends on A. Indeed, the 
state is not scale free - the wave function in the relative 
coordinate of the two fermions dies exponentially with a 
scale A^^ as noted in reference [211] . 

For a generic GFC, it is known that the critical scat- 
tering length required to induce a bound state is 

negative and finite^*^ and is given by Qsc = where 
is a dimensionless function. For a given < 0, 
this corresponds to a critical gauge coupling strength 
^(A) 



A. = 



The crossover with increasing A is then 



governed by the relative magnitudes of Xt and Ac. If 
Xc Xt, the crossover regime coincides with the regime 
of the FSTT. On the other hand if A^, ^ A^, the crossover 



regime is centered around A ~ Ac. In any case, for 
A >• max (At, Ac) the ground state will be a condensate 
of rashbons determined by the GFC in question. It is 
evident that except for the EP GFC, every other GFC 
will support a BCS-rashbon BEG crossover. 

We now discuss the situation with a small positive scat- 
tering length with kpUs <^ 1. In absence of a gauge field, 
the ground state is BEG of bosonic pairs of fermions with 
mass twice that of the fermion mass. In the presence of 
the gauge field, this BEG will evolve to the rashbon BEG 
as A — >■ 00, i.e, there is a BEC-rashbon BEG crossover. 

The authors are not aware of any experimental realiza- 
tion of synthetic gauge field in fermionic systems. The 
natural question that arises is if the parameter regime of 
A > At with a high symmetry GFC can be realized in 
experiments. We do hope that our paper provides the 
motivation for this direction of experimental research. 



Acknowledgement 

JV acknowledges support from GSIR, India via a JRF 
grant. SZ is supported by NSF DMR-0907366 and 
DARPA Nos. W911NF-07-1-0464 . VBS is grateful to 
DST, India (Ramanujan grant) and DAE, India (SRC 
grant) for generous support. We are grateful to Tin-Lun 
(Jason) Ho for discussions. 



* liayantha@physics.iisc.ern et.in| 
t shizhonK.zhanK@Kmail.corrT 

* |shenoy@physics.iisc.ernet.in 

^ Y.-J. Lin, R. L. Compton, A. R. Perry, W. D. P hillips, 
J^V. Porto, and I. B. Spielman, |Phys. Rev. Lett. 102, | 



130 401 (2009) 

^ Y.-J. Lin, R. L. Compton, K. Jimenez-Garcia, J. V. Porto, 
and I. B. Spielman, Nature 462, 628 (2009) 



Y.-J. Lin, K. Jimenez-Garcia, and I. B. Spielman, [Nature 



[~471, 83 (2011) 

K. Maeda, G. Baym, and T. Hatsuda, Phys. Rev. Lett. I 
[~10 3, 085301 (2009) 
° D. Jaksch and P. ZoUer, New Journal of Physics 5, 56 

(2003). 

® K. Ost erloh, M. Baig, L. Santos, P. ZoUer, and M. Lewen- 
stein, Phys. Rev. Lett. 95, 010403 (2005) 
J. Ruseckas, G. Juzelitinas, P. Ohberg, and M. Fleis- 

chhauer, Ph ys. Rev. Lett. 95, 010404 (2005); 

F. Gerbier and J. Dalibard, |New Journal of Physics 12, | 



033 007 (2010y, 

" See commentary on reference [3] by T.-L. Ho at 

http: / /www. condmatjournalclub.org. 
^° W. Ketterle and M. W. Zwierlein, "Making, prob- 
ing and understanding ultracold Fermi gases," (2008), 
,arXiv:0801. 2500 [cond-mat] 



I. Bloch, J. DaUbard, and W. Zwerger, Rev. Mod. Phys.| 
80, 885 (2008)| 

S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Phys. 80, 1215 (2008)!! 



Rev. Mod. 



A. J. Leggett, in Modern Trends in the Theory of Con- 
densed Matter, edited by A. Pekalski and R. Przystawa 
(Springer- Verlag, Berlin, 1980). 

A. J. Leggett, Quantum Liquids: Bose Condensation and 
Cooper Pairing tn Condensed-Matter Systems (Oxford 
University Press, 2006). 

^® P. Nozieres and S. Schmitt-Rink, Journal of Low Temper- 
ature Physics 59, 195 (1985). 

^"^ M. Randeria, in Bose-Einstein Condensation, edited by 
A. GrifBn, D. Snoke, and S. Stringari (Cambridge Uni- 
versity Press, 1995) Chap. 15. 

T.-L. Ho and S. Zhang, "Bose-Einstein Condensates in 
Non-abelian Gauge Fields," (2010), arXiv: 1007.0650, 
[cond-mat. quant-gas][ 

C. Wang, C. Gap, C -M. Jian, and H. Zhai, |Phys. Rev. 



I Lett. 105, 160403 (2010) 

J. P. Vyasanakere and V. B. Shenoy, |Phys. Rev. B 83, 
J)94515 (2011) 
E. Braaten, M. Kusunoki, and D. Zhang, jAnnals "of 
Physics 3"237 1770 (2008) 



The quadrupole operator 

1 (5-5/3 ^ gPgc.^ 



defined 
where S" 



D. M. Eagles, Phys. Rev. 186, 456 (1969) 



^^-^fJ"*^ where 5" are spin opera- 
tors. A nematic state has (5") = 0, while (Q"") / 0. 
In the present context, the singlet piece of the two body 
wave function does not contribute to the quadrupole mo- 
ment, while the triplet wave function has (5°) — 0, but 

{Q"") ^ 0. 

It is to be noted that the result for the EO case presented 
in [20] has a minor error of a factor of 2. The correct 



11 



results are quoted here; see also Section III B of e-print 
arXiv:1101.0411v2. 

Throughout the paper, the chemical potential is referred 
to the bottom of the + helicity band. 



R. B. Diener, R. Sensarma, and M. Randeria, Phys. Rev. 
A 77, 023626 (2008) 

J. P. Vyasanakere and V. B. Shenoy, (2011), under prepa- 
ration. 



